Preprint typeset in JHEP style - HYPER VERSION 



FTPI-MINN 07/02 
UMN-TH-2535/07 
MCTP-07-08 



SO(10) unified models and soft leptogenesis 



E. J. Chun 

Korea Institute for Advanced Study, 

207-43 Cheongryangri-dong, Dongdaemun-gu, 

Seoul 130-012, Korea 

e j chunOkias . re . kr 

and 

Michigan Center for Theoretical Physics 
Department of Physics, University of Michigan 
Ann Arbor MI 48109, USA 

L. Velasco-Sevilla 

W. I. F. Theoretical Physics Institute, 
University of Minnesota, 

116 Church St SE, Minneap olis, MN 55454, USA 
lilianaOphysics . umn . edu 



Abstract: Motivated by the fact that, in some realistic models combining SO(W) GUTs 
and flavour symmetries, it is not possible to achieve the required baryon asymmetry through 
the CP asymmetry generated in the decay of right-handed neutrinos, we take a fresh look 
on how deep this connection is in S'O(IO). The common characteristics of these models 
are that they use the see-saw with right-handed neutrinos, predict a normal hierarchy of 
masses for the neutrinos observed in oscillating experiments and in the basis where the 
right-handed Majorana mass is diagonal, the charged lepton mixings are tiny. In addition 
these models link the up-quark Yukawa matrix to the neutrino Yukawa matrix Y u with 
the special feature of Y{[ — > 0. Using this condition, we find that the required baryon 
asymmetry of the Universe can be explained by the soft leptogenesis using the soft B 
parameter of the second lightest right-handed neutrino whose mass turns out to be around 
10 8 GeV. It is pointed out that a natural way to do so is to use no-scale supergravity where 
the value of B ~ 1 GeV is set through gauge-loop corrections. 
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B.3 Results and comments 



1. Introduction 

Some of the unanswered questions left by the Standard Model (SM) acquire a new light 
in Grand Unified Theories (GUTs). Among them the understanding of the structure of 
masses has an important role. Theories based on 5*0(10) seem to have a special place 
due to the possible structure of masses that can be achieved. First of all the spinorial 16 
representation can accommodate all the known fermions including a right handed neutrino 
component, N, which makes it possible to embed in a natural way a see-saw mechanism for 
the explanation of the tiny masses of low energy neutrinos. On the other hand baryogenesis 
through leptogenesis Q is a simple mechanism to explain the observed baryon asymmetry 
in the universe. Here a lepton asymmetry is dynamically generated and then converted 
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into a baryon asymmetry due to B + L violating sphaleron interactions. The lightest of the 
right-handed neutrinos is produced by thermal scattering after inflation. It subsequently 
decays out-of-equilibrium to a lepton and a Higgs doublet producing a CP and lepton 
number violating asymmetry. The connection of leptogenesis to SO (10) GUTs is then 
natural, both by the inclusion of right handed neutrinos in the 16s and by the possible 
variation of their masses which can be some orders of mag nitude below M GU T ~ 10 16 GeV. 

However in realistic models where the expansion parameter describing the Yukawa 
couplings of neutrinos, e u is of the order of the expansion parameter describing the Yukawa 
couplings of up-type quarks, e u ||, || J|], the value needed for the mass of the lightest 
right-handed neutrino, Mjvi ~ lO^ 7-8 ) GeV, lies below the bounds for a successful thermal 
leptogenesis; Mjv x > 10 9 GeV ||. This prompts us to address two questions as follows: 
a) How general is this statement in the context of Grand Unified models where one naturally 
gets 



- r , . / : (1-1) 

VTL r 



b) If that is the case for many classes of such models, what are the possible scenarios for 
leptogenesis that we may require to consider in order to preserve such a feature of Grand 
Unified Models? 



In the present work we show how easy it is to generate the relation Eq. (1.1) and at 
the same time obtain Mn ~ lCr 7-8 ) GeV for lighter right-handed neutrinos. This will be 
shown explicitly in the limiting case of the vanishing neutrino mixing angle associated with 
reactor experiments s%3 — > 0, and of the dominant contribution of two light right-handed 
neutrinos. The contribution from non-zero S13 and the heaviest right-handed neutrino can 
be taken as perturbations. This enables us to determine, from the current experimental 
values, the possible form of the Dirac couplings of left-handed neutrinos in the basis where 
the mass of the right-handed neutrinos is diagonal, but without assuming any particular 
hierarchy among them. 

We then embed these results in a 50(10) context, without fully specifying a model 
but rather making choices that are compatible with GUTs, which can be used as a starting 
point in the construction of a complete SO (10) model. In fact our results are compatible 
with the models of g (RV), g (CM), § (BKOT), and §] (DHR), in which one considers 
symmetric matrices at the scale where SO(10) has not been broken and then explains 
the generation of fermion masses with a minimal content of Higgs bosons, such as two 10 
representations, which are the Higgs bosons in the u and d sectors, a 126 representation 
to generate masses for right-handed neutrinos and a non-renormalizable operator 45 to 
distinguish some features of the fermion masses. These are generic features that may be 
used to construct more specific models, and in fact have been used extensively (P]-Jl2|). 

The neutrino Yukawa texture and the right-handed neutrino masses determined in 
many works do not allow for the standard leptogenesis JTJ to produce the required baryon 
asymmetry of the Universe. Other mechanisms of inducing the appropriate baryon asym- 
metry in extensions of the MSSM may be implemented (e.g. |l3|] and see [jfj]] for a review) 
or also mechanisms just using extra right-handed neutrinos (e.g. p3[|). 
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However, we observe that the soft leptogenesis of (15], allowing a resonance condition 
with a small B term, works quite well within our scenario through the decay of the second 
lightest right-handed sneutrino whose mass is around 10 8 GeV. We will see that the required 
small B term (B ~ 1 GeV) arises from the gauge one-loop correction involving a heavy 
GUT gaugino once its tree-level value vanishes as it can be arranged easily in no-scale 
super gravity. 



2. Possible forms of vtl v lr and M R in models with an underlying Grand 
Unified theory 

We will identify elements of the right-handed neutrino mass matrix Mr and the Dirac neu- 
trino mass matrix, m u LR = Y v vsm (3/V2, which are related to the low energy observables 
(neutrino masses and mixing angles) by 



m 



U 1 



m 



m,. 



U* = -ml R M R 1 (ml R ) J 



(2.1) 



Here U is the neutrino mixing matrix and m Vi are the neutrino mass eigenvalues. This 
expression is valid in the basis where charged leptons are diagonal, if their matrix is not 
diagonal then we get U = U u U e * . Following the standard parameterization, let us write 



U = U 23 P^U 13 P 5 U 12 P m 



(2.2) 



where Uij is the rotation matrix in the plane, P$ = diag[e~ 4<5 , 1, e 4<5 ] and P m = 

diag[e v , e icr , 1] are the Dirac and Majorana phase matrices, respectively. We will use the 
notation {rn u LR )ij = niij and (Y v )ij = Uij, that is, mij = yijV sin/J/v^- 

For numerical values, we use the latest values of the fits for neutrino oscillation ob- 



servables [16]: 



Am^ G [7.3, 9.3] x 10 _5 eV 5 
Am| 2 e [1.6,3.6] x 10" 3 eV 2 



t? 2 G [0.28,0.60] 
4 G [0.25,2.1] 
s? 3 < 0.041. 



(2.3) 



Following the GUT relations consistent with the hierarchical pattern of all the fermion 
masses, we will work with the normal hierarchy; 



(2.4) 



Then, we make the following definitions to help us trace the hierarchy of low-energy neu- 
trino observables: 



m 



m 



(2.5) 
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2.1 GUT considerations 

Let us recall that the Yukawa sector of 50(10) at the renormalizable level comes from the 
following allowed couplings in the matter Lagrangian ]Jl7]] ; 

Cm = yjf (16)^16)^10) + 4 20 (16) i (16) i (120) + y}j 6 (16^(16),- (126) (2.6) 



due to the famous decomposition 16 © 16 = 10 s © 120 s © 126 s . As has been widely 
stated, the minimal Higgs content in order to generate Dirac masses for quarks and leptons 
and Majorana masses for right-handed neutrinos is to have two Higgs bosons in the 10 
representation of 50(10) and a Higgs in the 126 representation. 

The Yukawa couplings yj® associated with 10 representations give 

( m d ) T = m e , m u = m u LR , (2.7) 



while the 126 gives mass only to right handed neutrinos. 



Prior to the Super-Kamiokande experiments [18], there were successful models which 
could reproduce fermion masses with this Higgs content, however without explanation of 
the exact relations between the elements of the same Yukawa matrix. With the current 
data of the neutrino oscillation experiments [flq] , now it is clear that the minimal content 
must be extended to fit neutrino masses. This can be done by adding more 10, 120 or 
126 Higgs fields or by adding non-renormalizable operators or other elements in the theory 
beyond a GUT. 

If we consider non-renormalizable operators we can alter the structure of the Yukawa 
matrices. The Lagrangian of these operators is 

CM(ht)=O ii V°i, 0u = 16;#---#io|^...#16 j , (2.8) 

M[M) ytr 3 Mi M k M k+ i M e J ' y J 

where i?'s are possible representations coupling to 10 and 16 and yf- is their corresponding 
Yukawa matrix. An interesting case is the adjoint representation 45 whose vacuum expec- 
tation value (vev) can point to any direction in the space spanned by the 45 generators 
as long as it leaves the SM group SU(3) C x SU(2)l X U(l)y unbroken. The space of 45 
vev, that leaves unbroken the SM group, lies in the two dimensional subspace of U(l)' 
generators of 5*0(10) that commute with 5lV(3) c x SU(2)l x U{1)y- For example when 
5O(10) is broken to the SU(3) C x SU(2)r x SU(2)l subgroup, the general form to which 
the 45 should be pointing is 

(45) = (B-L + kT m )M 45 . (2.9) 

Successful and predictive scenarios can be obtained when incorporating flavour symmetries 
into GUTs. With a clever choice of a flavour symmetry, it may not be needed to invoke 
additional Higgs fields beyond two 10s and one 126 representations. However majority of 



such examples do require at least one non-renormalizable operator (see [19] for a review). 

From the experimental information on Vqkm an d the quark masses we can reconstruct 
the possible forms that the quark mass matrices, m q , can acquire from the above consid- 
erations. However since we are just able to construct the left diagonalizations of m q from 
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Vckm w e need to make more assumptions on m q . One of the most successful and well 
established assumptions is to consider that the mixings in each quark sector are small such 
that they remain small when combined in the Vckm- This assumption can be naturally 
identified with a strong hierarchical structure of m q . 

The diagonalizing matrix of Y q = \[2m q jv q can be parameterized as the multiplica- 
tion of unitary matrices U 23 Ui 3 Ui2 where each Uij is a rotation matrix in the sector 
including some phases. When Y q is hierarchical, the angles corresponding to this pa- 
rameterization can be in fact approximately identified with the angles obtained from the 



approximate multiplication [2' 
the up quarks, we have 



U1 2R U1 3R U1 3R YU 23L U 13L U 12L 



lag- 



For the case of 



Du ~ I c 12 c 122/ll 



„i<t>L r L R II 



e *^' Rc 12' s 122/21 + ,s 12' s 122/22l 



L R II , i<j>L r L R 11 1 p i<t>R\„R „L II , n u ,, 1 
c 12 c 12tf22 ^ e L 12*122/l2 t e I L 12 6 12y21 ^ 6 12 6 12i/ll I 



~.R „L 



Vt « I2/33C23C23 



(2.10) 



where y" is the matrix after the unitary transformation in the (2,3) and (1,3) sectors. 
We are particularly interested in analyzing the behavior in the (1,2) sectors. For small 
rotations we have cf 2 = c§ k, 1, sf 2 = y 2 i/y 22 and s{ 2 = y'{ 2 /y 22 . When y'{ x < 2/122/21/2/22; 
one gets 



\yu\ 



I s 12 s 12S/22 1 



\JJc 



I2/22I 



(2.11) 



which, for the symmetric case, gives the famous relation sf 2 
basis of the Gatto-Sartori-Tonin relation f^2]l 



yj m u /m c . This forms the 



V,, 



d i$i u 

6 12 e 6 12 




(2.12) 



This relation is in good agreement (originally the contribution just from the down sector 
was assumed, but now the contribution from the up sector has become important) with the 
experimental value of V us . In order to review the conditions for which y xl <C y'\ 2 y 2 \l y 22 -, 
let us give explicitly their expressions in terms of the small rotation angles and the original 
Yukawa couplings: 



yli 



y'k 



yu 



4, 



+ 4yi3][s232/2ie^23 + c % 3 y 31 ] 



1 oh 



y33C2 3 cf 3 + 0(2/232/32/2/33) 



[c|2/12 



sl\y 13 e *^][<4lft2i 



S232/3ie 



c 23 s 232/32e" 



?4 



c 23 S 232/23e + 5235^37/33 



-10, 



(2.13) 



From these expressions it is clear that, if the original Yukawa matrix contains a zero in the 
(1,1) position, then the inequality 2/11 <C 2/122/21/2/22 is immediately achieved because of 
2/33 ^ 2/22 leading to the desired hierarchy between charm and top quark masses. However 
notice that 2/11 does not have to be exactly zero. As long as it is suppressed enough with 
respect to 2/122/21/2/22 we can always have a relation of the type Eq. (|2.12| ). This condition 
is referred hereafter as the limit — * 0. As we have seen it can be realized in a particular 
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basis of the Yukawa couplings satisfying Eq. ( |2,11| ), which is however a basis independent 
statement. 

For a symmetric case, if we assume yx2 2/23 then in any of the cases of having 
2/22 ~ 2/23, 2/22 2/23 or J/22 > 2/23, the second term in the first of the equations of Eq. ( 2.13| ) 



will be smaller than the expression for y'^- Thus the interesting case is when yn is the 
leading contribution in y'[ x . Let us take the case j/22 ~ 2/23 then the leading contribution in 
2/12/2/22 W1U ^ e simply 2/12/2/22 and hence we obey the condition y'[ x <C 2/122/21/2/22' m terms 
of the original matrix elements: y\\ <C 2/12/2/22- 



For matrices of the form [20], |2J, |23|], where the elements Y^ and Y 2 2 are respectively 



and e„, one simply needs to require Y^ <C e„. For matrices of the form 24], which 
are also symmetric and hence compatible with SO (10), one recovers the requirements of 
the element Y^ with the present analysis. In this texture, Y^ 2 is e„ ~ A 6 , Y^ is zero but 
Y23 is e„, and hence (y") 12/2/22 = e4 - Thus, by making Y x " = 0(e„) in this case, an 0(1) 
correction to the relation s" 2 ~ \l Tfhjrnc can be made. 

2.2 Compatibility with the experimental information 

Symmetric fits to the quark masses can be used as a guideline to construct models with 
underlying 50(10), or SU(4) C x SU(2) L x SU(2) R , which after the 5O(10) breaking assume 
a symmetric structure. However a small departure from symmetric matrices do not change 
the qualitative behavior of such fits and can be made compatible with 5O(10) and flavour 
symmetries. 

The fits of these matrices into the experimental information can be made in many 
ways depending on our theoretical assumptions. A minimal choice is to assume that the 
supersymmetric corrections to the quark masses will not have a strong impact on the ratio 
of masses that we use for the fit. We also assume that the Yukawa matrices are the only 
source of CP violation and that the contributions from the transformation of the squared 
soft mass matrices of the supersymmetric particles to the basis where the Yukawa matrices 
are real and diagonal are negligible. We specify this last requirement because we can choose 
to fit the quark mass matrices with ratios of masses and with the fits of the unitary triangle, 
where various CP violating experiments are taken into account. 

With increasing precision in the determination of the fits of the unitary triangle in the 
SM, however, one has a very tight constraint on the parameters and one should not regard 
this as a final fit of a particular texture but as an indication of the current compatibility of 
such texture with the theoretical assumptions and experimental information. The purpose 
of the present analysis is to clarify the consequences of having a negligible Y^ entry, giving 
the relation s" 2 ~ \/ m u/ m c and furthermore a low range of right handed neutrino masses, 
which turn out to be incompatible with the standard thermal leptogenesis to produce the 
observed baryon asymmetry of the universe. For completeness we present in Appendix ([B]) 
the current fit of a symmetric texture with negligible Y^ element, compatible at 68% C.L. 
with up-to-date fits of the unitary triangle. 

It is interesting to see that two different choices of symmetric matrices in the up sector, 
((2|, U, j|], U) and (]24f|, p5|), give rise to different phenomenology and consequences, 
e.g., for leptogenesis. This is because we then have clearer selection criteria on how to single 
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out the models, apart from direct fittings to experimental mass matrices and quark mixing. 
Further precision analyses in the unitary triangle will constrain more such possibilities. 



For non-symmetric fits [21], we have a rather different situation because often tb'q 
should be of the order of m\ 2 to reproduce appropriate mixings and quark masses. The 
non symmetric fits can be used to construct models with underlying SU(5) symmetries 
(see for example p6|]). 

In the following, we will consider the limit (w-^)n = nifi — ► and explore its conse- 
quences for neutrino oscillations and leptogenesis. 



2.3 Consequences of (m1 R )ii — ► 

2.3.1 General structure of the matrix of low energy neutrinos 

The explicit form of the low energy neutrino mass components from Eq. (|2.1|) is 



m 
m 



V3 
V 

efj, 



m 
m 



V 

er 



m 
m 



V3 
V 



in 



/XT 



rn 



V3 



to; 



in 



2 2 , U{8-a) 2 , t 2 2 -2ip 
lc 13 6 12^ e 6 13 + lc 13 c 19 e 



13 



-13 12 t 



S23Sl3Cl3e l(l5 2<J) + rCi3Si 2 (ci2C 2 3 



512^13^23 - 
c i5 



^-C 12 C 13 {C 2 3S 12 + — 5 ) 



i(S-2a) /C23S12S13 . s . t C12S23S13, 
C23Cl3Sl3e V ; - rc 13 S 12 { — h h C12S23) + ^Ci 2 Ci 3 (C23Si2 ~ h ) 



2 2 -2ia 
s 23 c 13 e 



+ r(ci2C 23 



Sl2Sl3S23 N 2 



2ip 



(C23«12 + Cl 2 Sl3S23^ 



C23S23C13 ,C 2 3Sl2Sl3 . w Sl2Sl3«23 N 
r( ~ h Ci2S23j(Cl2C23 ~r ) 



2io 



AS 



t , , Ci 2 S23Sl3 N , 
^(C23S12 + — s )(S12S23 



C12C23S13 
c i<5 



cf 3 c| 3 ,C23Sl2Si3 2 t 

+ r ( —S + C12S23) +^ 2 -(s 12 S 2 3 



a 2ia 



ci 2 c 23 s 13 2 

e i8 I ' 



(2.14) 



where we have expressed the elements of m v in terms of the sub- indices e, [i and r. In 
Appendix (Eh we have written the numerical central values of the angles, up to t, sf 3 and 



possible phase variations. From Eq. (A..1) we can see that in the limit of S13 = 0, all 
the numerical entries that are not multiplied by t in Eq. ( [2.14 ) do not change its order 
of magnitude. On the other hand the complete form of m u in terms of a diagonal matrix 
Mfl = diag[Mi, M 2 , M3] and a general matrix m v LR is given by 



m 



y- 

4- Mi 



in 



Li 



™li™2* m li m 3i 



<2i 



m 



2/ 



to 9 ,-to 



2i" l 3i 



m \i m 3i m 2i m 3i m 3i 



(2.15) 



Now when {rn v LR )n = mu — ► 0, this matrix acquires a very simple form. Then we can 
simply identify the elements of Eq. ( 2.14 ) with Eq. ( p. 15 ) and find the form of Mj in terms 



of rriij and the restrictions of its elements. 



It is clear from Eq. (|2.14| ) that the contribution from t can become relevant just for 
But since t < O(0. 1) according to Eq. ( |2.3| ), this contribution can be at most 



mg e and m v eT 
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as the same order of the rest of the contributions in (m u )n or {m v )\^. From Eq. ( 2.15 ) we 
can see that these elements are given by 



m\ 2 



+ 



771 



13 



mi 2 m 32 mi 3 m 33 



(2.16) 



M 2 M 3 ' eT M 2 M 3 

Identifying certain elements in m v LR , we can obtain the predicted ranges for the right- 
handed neutrino masses which will be derived in detail in the next subsection. In this 
subsection, we try to extract the general expressions for the parameters such as 

M 1 Mi _ (Y^Y u ) n v 2 

or arm m.i = 

M 2 



Mi, 



Mi 



Mi 



(2.17) 



which are relevant for leptogenesis. Note that m" e is the element which contains less 
parameters and so we can make less assumptions when deriving expressions for Eq. ( 2.17] ) . 
Then we have 



M 2 



2 ri , Mi-rn^ 1 
m 12L i + M 3 mf\ 



,[rc? 3 s? 2 + e»(*-^) a 2 3 + tcf 3 cf 2 e- 2 ^] 



■/n, 



> 2 x 10 16 y 2 2 



sin 2 /3 



1 + 



M 2 mf 3 
M 3 m 2 



12 



GeV, 



(2.18) 



where we have set the bound on M 2 by taking the numerical values of the first two contri- 
butions of the denominator in Eq. ( 2.18|) . 

From this relation we can study the behavior of the parameters in terms of 



M 2 m 2 3 
M 3 m 2 



(2.19) 



12 



For a <C 1 we can see that the order of M 2 is determined just by y\ 2 . In this case, 
considering the expressions of mj^ + m^ T , we obtain the ratio of Mi to M 2 given by 



Mi 
M2 



m u ee {m 2 2 i +m§i) 



(r< +r< T )m 2 2 (l + a) 



a(m 2 3 + m 2 3 ) 



t 2 - + (m\ + m 
*i3 v 61 



(2.20) 



Then, we find that generic 50(10) models lead to rh\ given by 

mi = max|m7,„ + m^, b m'j , where 



m ax[m^ + m v TT ,b m v ee ] , 



b = 



a(ml 3 + m 2 3 )^ 



+ (m 2 2 + m 



22/ 



mf 2 (l + a) 



(2.21) 



We can consider the cases for b < 1 and 5 > 1. For 6 < 1, the order of magnitude of rhi 
is fixed simply by + m rr ~ (10~ 2 , 10" 1 ) eV. For b > 1, rhi ~ 6(10~ 3 , 10~ 2 ) eV. This 
points out an important consequence that the right-handed neutrino gets out of thermal 
equilibrium when it is very non-relativistic [see next section] because < rhi where 



16vr 5 / 2 i /2 v 2 

m * = 7^9* TT- 

3^ M P 



1.6 x 10 - - 3 eV 



(2.22) 
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with = 225 for the relativistic degrees of freedom in supersymmetric standard model. 
This brings a strong suppression to the CP asymmetry generated through the decay of the 
right-handed neutrino. 

For a>l, instead of the bound for M2, we determine the bound for M3, which can be 



obtained by taking the replacements: M2 <-> M3 and 77143 m i2 m Eq. ( 2.1 8| ) . Similarly, 
the corresponding rh\ is given by the relation ( |2.21| ) with exchanging the indices 2 «-> 3 for 
the expression of b. 



This has the same behavior as Eq. ( 2.21| ) except that in GUT models the natural 



assumption is to have 2/33 = 0(1) and hence the second factor bm v ee is likely to be the dom- 
inant. Then this contribution to rh\ goes like rh\ ~ m ee y^/y^ which could be significantly 
bigger than m* ~ 10~ 3 eV. 

In all the models p|, p[, H and |6|] in the basis where Mr is diagonal 1 the condition 



Eq. ( 2.11 ) is satisfied for the case of the Dirac coupling of neutrinos and then we can 
identify its behavior in terms of a and b. In [||, a particular realization of the case with 
a 1 and 2/33 = 1 was explored and one gets Mi ~ 10 7 GeV and the wash-out factor, 
~ 10 -5 , to the leptonic CP asymmetry in the decay of right-handed neutrinos. The other 
parameters can bee seen in Table |. In §], one gets a = because effectively only two 
right-handed neutrinos are taken into account and the mass of the lightest right handed 
neutrino is Mi = O(10 7 ) GeV. In |, a « 1 and b G (O(0.1),l) and the mass of the 
lightest neutrino is Mi = O(10^ 7,8 ^) GeV. A realization of the case with a <C 1 and b ~ 1 
was explored in |7j and || where it is also not possible to achieve a successful thermal 
leptogenesis, due to the washout factor, although the mass of the lightest right handed 
neutrino is of order Ml = O(10 10 ) GeV. In Table [l] we have summarized the properties 
relevant for leptogenesis of these models and we will make more comments about them in 



Section 2.3.3 



In the next section we take the limit 

si 3 ^0, t->0 and 1/M 3 -»■ 0, (2.23) 

which gives a <C 1 and b <C 1 and allows to understand more deeply the connection of 
( m Li?) 11 ~^ with the low mass of Mi in GUT models. 

2.3.2 Limit of S13 — > 0, t — > and contributions proportional to I/M3 negligible 

The goal of this analysis is to identify the possible shapes of m v LR and Mr in the basis 
where Mr is diagonal and the form of the mixing angles in terms of their entries. We will 
assume that the mixing of charged leptons is small and hence can be ignored. In the limit 
under consideration, Eq. ( |2.14| ) becomes much simpler as follows: 



S?2 r C12C23S12I" -Cl2«12S23r 

ci2C 2 3Si2r cf 2 c| 3 r + sl 3 e~ 2ta -cf 2 C23«23r + c 2 3S 2 3e~ 2la 

-Cl2Sl2.S23r -Ci 2 C 2 3S23r + C23.S23C" 2 " 7 cf 2 S23 r + e ~ 2l ° C 23 



. (2.24) 



1 For the models pi, |^|, this transformation has been performed since for those models Mr is not diagonal 
in the basis where the underlying symmetry is broken. 
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Model Ref. 



M fl [GeV] 



M 2 mf 3 
M3 m'i 9 



mi 



BKOT 



CM 



DHR 



RV 



10(7,8) 



10 9 



io( 



14,15) 



10 7 



10 a 



10 



12 



10 



10 7 



10 



10 1 



10 



13 



10 8 



> 10 



10 



II. ~ 10" 3 

III. ~ xl(T 4 3 x 10 2 

IV. ~ xlO" 4 



10 



10" 



10 3 

O(l) 

O(l) 



b m u ee ~ 
(10 _1 ,1) eV 

b m u ee ~ 
(1, 10) eV 



m U3 < 0.625 eV 



Table 1: Models based on ,50(10) which do not generate the observed amount of baryon asymmetry 
through the decay of the right handed neutrinos in thermal leptogenesis. 



When comparing the elements m^- of Eq. ( p. 24 ) with those of Eq. ( p. 15 ), we have six 
equations to solve, which are more than the low-energy parameters to determine: one mass 
ratio, two angles and one phase. Thus the elements of mP and Mr are more restricted. 
By comparing the m v ee component in Eq. ( |2.24j ) and Eq. ( |2.15| ) , we can readily identify M2 : 



Mo 



m 



12 



m„„ rs 



Considering the ratio m u eT /m y e ^ we obtain an important relation: 

m.32 



*23 



rn-22 



(2.25) 



(2.26) 



Analogously the ratio m^ r /m^ r leads to 
*23 = 



m 21 m 3 l + fm 22 m 32 



m 31 + M2 m 32 



1 + c 2 2 t 2 3 re 2 ^ 
1 - cf 2 re 2i<J 



(2.27) 



When we determine the ratio Mi/M 2 we can put a restriction on the elements 77721 and 
771,31 from Eq. ( |2.26| ) and Eq. ( |2.27| ). Now the solar mixing angle is given by the following 
simple equation; 



*12 



77712 



"732 



"^22 \Jm\ 2 + 777 



(2.28) 



:>>2 



which is obtained by considering va v ee jvn v etl and the relation i 23 = m 3 2/ m 22- Adding 777 



and m v TT we obtain 



Mi 
W 2 



rs 



2 1 

12 



m 



12 



m 21 + m 31 



+ rc 2 2 (l-i 2 3 ) 

















.™12. 



{m\ x + 777 



2iu 



31) 



(2.29) 
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where the last equality follows from Eq. ( 2.28j ). Analogously by adding 
obtain 



m^ T and m v TT we 



Mi 




" r 4T 
m i2. 


[m 21 m 31 + ml] 


M~ 2 ~ 


e~ 2ia c 23 (c 23 + s 23 ) + rcf 2 | 


S 23 (s 23 - C23) - ^[ m 32 + "l22m 32 ]) 



(2.30) 



Now dividing Eq. (|2T29|) by Eq. ( p3o| ) we can find solutions for p = m 2 \ / 71131 and then all 
the parameters can be expressed in terms of two unknowns; 77131 and m 22 . It is illustrative 
to consider the limiting case of i 2 3 = 1 given the fact that the atmospheric neutrino mixing 
is the best measured quantity; t| 3 = 17 ' 3 - Note in this case that we have the two solutions 
p = or p = 1 along with the following relations; 

m 32 = -m 22 and m\ 2 

Equating Eq. ( gjg ) and Eq. ( gj27fl we find 



2ti 2 m 22 • 



(2.31) 



Mi 
M2 



771. 



31 



1 



7/1 



fc with k =-[(l-p)-(l+p)ci 2 re^], 



22 



which leads us from Eq. ( p. 25 ) to 



Mi 



-2m 



31 



fc. 



12 



rm,, 



(2.32) 



(2.33) 



Now we find that Eq. ( |2.32D is indeed compatible with Eq. ( 2.2S| ) for p = 1. 



Form of the matrix m£ R 

Allowing a deviation from the limiting case of t 23 = 1, we can get a more general values 
for the parameters. In any case, we note that given m 22 we can determine the elements 
77i3 2 and 77ti2 through low energy observables; t\ 2 and 7j 23 as in Eq. ( p. 28 ), and hence fix 
the scale of M 2 by using m U3 and r. Although we cannot fix the values of m 2 i and m 3 i 
independently, we can fix their ratio p and satisfy all experimental constraints. Then we 
can see how the hierarchy of Mi and M 2 depends on the ratio of and 711,31 . 

That is, we determine the following structure of m v LR 

tl2, 



m LR 



-^77122 

C23 ZZ 



pm 31 77122 x 2 , (2.34) 
"131 ~t 23 m 22 x 3 

where the elements of the third column cannot be determined due to the limit I/M3 — > 
we have taken, and we can see that all the entries in the second column are of comparable 
order. As we mentioned before, the motivation for having (711^)11 — > was linked more 
closely to the symmetric fits, which implies mi 2 ~ ?n 2 i. From Eq. ( |2.29| ), we can write 



Mi 
M2 



777. 



m , 



21 2 1 +P 
— rsi2 ^ 



(2.35) 



'12 P 

which leads us to a conclusion of M\jM 2 ~ 0.2 with p = 1. In the following, we will give 
more precise determination of the right-handed neutrino masses further considering the 
50(10) structure of m u LR in Eq. 
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2.3.3 SO(10) and the scale of Mi 

We started our discussion by considering an underlying SO (10) framework with the mass 
matrix relation (2.7). Then, it has been realized Q that at least one non-renormalizable 
operator should be included in the Yukawa sector to get a successful fit to the quark and 
charged- lepton sectors. Generically the representations which break the 50(10) down to 
SU{5) or SU(4) C x SU(2) R x SU(2) L and then further down to SU(S)c x SU(2) L x U(1) Y 
do not give rise to the appropriate structures of fermion masses, and thus we need to 
add at least one non-renormalizable operator. As we mentioned, the non-renormalizable 
operator involving the 45 representation is quite useful since its vev can be aligned in the 
B — L + /cT#3 direction for k £ Z whose value is different for every species of fermions. 
As we discussed in Section 2A, fits to the quark masses [pQl -|^3||, give a structure of the 
up-sector of the form 



U V w 

v x y 
w y z 



, u <C v, u < w < v, v < x, y < x <C z ~ O(l), (2.36) 



where the exact symmetry of the matrix is not necessary but the orders of magnitude in 



the elements Y£j and Y^ have to be the same. In the previous Section |2.3,2| , we determined 



the possible form of m u LR which can be naturally understood in the context of the SO (10) 
models when 77131 ~ m 22 • There are two forms of Yukawa matrices for the Dirac neutrinos 
that have been exploited in the literature: 



u 



V 



12 U»13 
v '2l x 2/23 
W 31 y'32 z ' 



v! <C w' 13 



- w' 31 ~ v' 12 ~ v' 21 , v' 12 ~ x' ~ y' 32 ~y' 23 <.l 

' W>31 < v 'l2 ~ u 21i v 'l2 < x ' ~ y'3,2 < y'23 < 1 



2/32 ~ ^23 < 1 



(2.37) 



Each of these has been justified within the context of S'O(IO) and a particular flavour 
symmetry. The models g} (BKOT), @ (CM) and |] (RV) are examples of the first form 
of Y u in Eq. ( ^.37|) . In for instance, the behavior of x 1 ~ w 1 ~ v', which is different 
from v <C y < x in the up-quark sector, was explained by the coupling of the elements (2, 2) 
and (2,3) to a 45 representation whose vev is proportional to B — L + 2Tr3- This vev 
is different from zero for up-quarks while it vanishes for neutrinos and thus one is forced 
to take into account the next leading contribution which has to be of the same order for 
the entries (1,2) and (1,3) leading to a large mixing solution. The model of § (DHR) is 
an example of the second case. It has also a coupling of a 45 representation in the (2, 3) 
sectors. But the difference with respect to the up-quark Yukawa couplings is given by the 
vev breaking the flavour symmetry D 3 x U(l) x Z 2 x Z 3 and together with the other orders 
of magnitude in Y v a large mixing is achieved. Now, Eq. ( p.!2| ) fixes the order of € u as 
follows; 



v 

X 



rrini 

— ~ (3,6 x 10" 



x 
z 



(2.38) 
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For the first case of Eq. ( fT37D , v' is straightforwardly related to v: 



v ~ e„ 



(2.39) 



On the other hand, for the second case of Eq. ( fr37D , u' is also related to the up-sector but 
due to the structure of D3 x {7(1) x Z2 x Z3 it is given by 

t/ ~ e 2 u . (2.40) 

Once the scale of mf 2 is fixed, we can determine the order of magnitude of M2 with = 
c(5 x 10~ 2 ) 3 , 



Mo 



9 9-9 

yf 2 ^ sin 



(3/2 



in. 



c 2 sin 2 /3 (0.5, 2.4) x 10 8 GeV, 



(2.41) 



where c is an O(l) number and va v ee is taken from Eq. ( |A.1| ) considering the normal hier- 
archical spectrum of low energy neutrinos, t <C r, and the uncertainties in Eq. ( |2.3P are 
taken into account. When r ~ t, non-trivial phases can make vn v ee very small and the above 
estimation of M2 has to be considered just as a lower bound. For Eq. ( p. 41 ) we have a <C 1 
and b <C 1, and hence this estimation can be applied to the BKOT, CM and RM models. 
For the model DHR || having yu = 6.27 x 10~ 3 = 0(e 2 ), one gets much larger value: 
M2 ~ (0.3, 1.1) x 10 12 GeV which agrees with the value quoted in Table |]. For the explicit 
case with p = 1 analyzed in Section 2.3.2 , the ratio of M\jM% becomes 



w 2 



in 



rii = ^|i(o.089,0.19), 



111 



(2.42) 



22 



Thus, we get the ranges of M x » (0.0045, 0.46) x 10 8 GeV. 



3. Baryogenesis through Leptogenesis 
3.1 Thermal Leptogenesis 

The CP violating asymmetry generated in the decay of heavy right handed neutrino into 
a Higgs boson and a left-handed lepton and its CP conjugated channel, and its supersym- 
metric counterpart are 

^N it i - r N j - Tftj 

1 N it l + 1 N t ,l L N l ,l + L N l ,I 

where r N i: i = £ aij gr(./Vj — ► l a H^) and T N _ j = Yl^pT^N^ — ► l a H^) are the Ni decay 
rates into I and I respectively. For the right-handed sneutrino JVj decay rates, one has the 
final states with the lepton I and slepton I. The B — L asymmetry generated by the right 
handed (s)neutrino decays is given by 

79 

Yb-l = -THfaYfi + e$Y^} = -Y B (3.2) 
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where Yg is the resulting baryon asymmetry converted from the Yb-l by the electroweak 
sphaleron processes, and rji is the efficiency factor that measures the number density of 
Ni/Ni decays at low temperature T « M ; and Y e i = Y eq (T > M;) = 135C(3)/(4tt 2 s*). 
For g± = 225 in the supersymmetric model, one gets Y^ q = 1.9 x 10 -3 . Since iVj and iVj give 
the same contributions in the supersymmetric limit; = = Tj and = e^. = e*, 
one can express the final baryon asymmetry as 

Yb = 1.3 x 10~ 3 rji 6i (3.3) 

where the observation requires Yb ~ 10~ 10 . We recall that the decay rate for the right- 
handed (s)neutrino is 

_ {Y^Y v D ) u M Ni 
Li-L Ni> l + L Niil - — . (3.4) 



Then the CP asymmetry can be expressed by 

(3.5) 



-y Im 




where, whenever the hierarchy Mir./Mjjr. = x is good enough, the function f(x) is just 

In our case, the right-handed neutrinos are charged under the SO (10) gauge group 
and thus are in thermal equilibrium at high temperature. The efficiency factor rji can be 
calculated given the value of Ki = Ti/H(T = Mi); 

K l = ^ = ™\ T (3.6) 

m* 1.6 x 10~ 3 eV v ; 

where the effective neutrino mass is defined by rhi = AirTiV 2 /Mf . If K{ < 1, the effi- 
ciency reached its maximum rji = 1. However, when Ki 3> 1 as is the case in most SO (10) 
models, the inverse decay remains effective for T < Mi and its decoupling temperature 
Zi = Mi/Ti is given by [§8| 



^zfe- z ^l + ^z t ^z t -l (3.7) 

and the corresponding efficiency factor can be well approximately by the simple form 

m ^^r(l- e~ z ^ /2 ) ■ (3.8) 
ZiKi V / 

It is important for our case to notice that the lepton asymmetry along the electron 
direction generated by the second lightest right-handed neutrino N2 is not washed out by 
the lightest right-handed neutrino Aq as we have yn —* 0. Therefore, let us consider the 
possibility of a successful leptogenesis either from Aq or N 2 whose effective neutrino masses 
are 

v 2 

fhx = — sin/3 2 |y 3 i| 2 [l + p 2 } rs m U3 « 0.05eV , 



rh 2 = sin f3 2 1 y 22 \ 2 1 + ^ 12 sa m V2 sa 0.009eV . (3.9) 
1W2 



-23 
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Thus we get 



<^)h£S till <"»> 



leading to 

' 8.5 x 1(T 3 for i = 1 



(3.H) 

6.4 x 10" 2 for i = 2 



On the other hand, one can readily check that we have 

3 3 M > - 



~ w*w s 10 <3 ' 12) 



putting the numerical values determined by the procedure of Section 2.2. This is too small 
to produce the required amount of baryon asymmetry, 10~ 10 , as expected from a general 
discussion with hierarchical neutrino mass spectrum [||]. 

3.2 Soft Leptogenesis 

It has been pointed out [15] that the 5-term soft supersymmetry breaking of the right- 
handed sneutrino provides an additional source of lepton number and CP violation, where 
the relevant couplings are given by 

-£ soft = {m\)iN%N Rj + ((ji e i Lj e* m h d + a^ N l Lj N* m h u + h.c) 

+ Q(M^Mi+h-c-)- (3-13) 

The effects of &jv = BMr terms are usually ignored because they are assumed to be highly 
suppressed by the difference in scales of the typical supersymmetric masses, 10 3 GeV, with 
respect to the masses of the singlet neutrinos, Mr > 10 7 GeV. It turns out that there is a 
region in the parameter space of B and Mr compatible with models for which the masses 
of right-handed neutrinos are as low as Mr ~ (10 7 — 10 8 ) GeV. The non-vanishing value 
of the generated lepton asymmetry is a pure thermal effect since at T = the generated 
lepton asymmetry in leptons cancels the one in sleptons: 

e- - ~-e- - 4 ^ [lmA] (314) 

i 

where A = a^/Y u . At finite temperature T / 0, the difference between the fermion 
and boson statistics yields non- vanishing lepton and CP asymmetry of the form €i{z) = 
e fj ^iH d &BF{z) where 5bf(z) can be approximated for z 3> 1 by the analytic function of 
5bf(z) = 2\f2K\(y/2z) / K\(z) with K\(z) is the modified Bessel function of the first kind 
@. 

Thus, in the soft leptogenesis, one gets the reduced efficiency defined by 



2^M^£> x m (3.15) 
Ki{Zi) 



- 15 - 



where z% is the decoupling temperature of the inverse decay calculated before. For our case 
with Zi = 7.51 and 4.96, we get 5bf{zi) = 0.105 and 0.3 and thus 

f 8.9 x 10- 4 for i = 1 ,„ 1c v 

L 1.9 x 10 ^ for i = 2 

Therefore, we require £f[._j[jj ~ 8.6 x 10 -5 and 4.0 x 10~ 6 correspondingly for a successful 
leptogenesis. When Ti < B, we get 

ImL4] _ rhjMj Im[A] 

e N^W d ~ B M . ~ 4^2 B ■ V-U) 

Now one can find that the lightest right-handed sneutrino N\ cannot produce enough lepton 
asymmetry due to a strong wash-out suppression. However, in the case of N% with M% = 10 8 
GeV, the desired value of tf l2 _ > i Hd ~ 4.0 x 10 -6 is found to be achieved for the hierarchical 
choice of soft parameters; Im[A] ~ 1.7 TeV and B = 1 GeV. Note that the leptogenesis 
scale ~ 10 8 GeV can be marginally allowed in view of abundant unstable gravitinos which 
decay late and upset the standard prediction of the big-bang neucleosynthesis . Recent 
analyses showed that the upper bound on the reheat temperature is Tr = 2 x 10 6 — 3 x 10 8 
GeV for the typical gravitino mass range of 7713/2 = 10 2 — 10 3 GeV, assuming the hadronic 
branching ratio of the gravitino decay is 10 -3 . We also remark that the bound on the reheat 
temperature can be loosened if the gravitino is stable and forms dark matter. In this case, 
the next lightest supersymmetric particle needs to be a stau and the reheat temperature 
up to Tr = 10 10 GeV can be acceptable plj . 



4. Supergravity description of a small B term 

It has been pointed out that the smallness of the B term may arise if we have for example 
a dynamical mechanism that sets B = at the leading order by arranging a specific form 



of the superpotential Ref. 32]. However, we find it difficult to achieve such a mechanism 
without introducing a fine tuning of parameters in the general supergravity context. In 
this scenario, the other ingredient to produce a small B term of the order m^^/^Ni it 



may be through a term / d6 A X^XN 1 Ni as in §§-§4 

To illustrate difficulties in a dynamical set up of B = in supergravity, let us consider 
the superpotential suggested in Ref. |32| 



W = n($i)NN + Af(&i)X + W', (4.1) 

where $j can be an observable field such as an multiplet of SO(10) (that is, 126 as iV can 
be 16) or any field in the hidden sector. The minimization condition of the scalar potential 
V = e K \K fi FiF-j - 3|H/| 2 ] reads Vi = where 



V, = e K 



KjfFiFj + [Fi(WK 3l ) + F 3 (W U + KtWi + K a W)] - SW t W + K t V (4.2 



is the derivative of the potential V with respect to a field /. Here we have defined Fi 
Wi + KiW. The b term coming from the scalar potential is given by 



e K 



K ij F 1 {( M + K if j,)-3(j,W + 2piW (4.3) 
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where the last term comes from \F X \ 2 = \2nN + Kj^W\ 2 with the minimal kinetic term for 
iV; Kn = N. On the other hand the minimization condition for X is 

= K*f [Fi(WK ]x ) + Fj(W iX + K t W x + K iX W)} - 3W X W, (4.4) 

assuming K x = and V = at the minimum. The indices i and j in the above equation 
contain $j and X. Assuming there is no mixing term between them in K (that is, K$. x = 
and K^ iX = 0), we separate the X index to write 

= K*J [F 3 {W iX + K.Wx)} - 3W X W 
+K xX [F X WK XX + F x (W xx + K X W X + K XX W)] . (4.5) 



Thus we need to arrange the second line of Eq. Q4.5|) and the Wix contribution to sum up 
to 2aW to cancel the above b term in Eq. ( ^4.3[ ). We find that, with an specific form of a 
Kahler potential, one can achieve such a condition which however requires a fine tuning of 
the parameters involved and does not have a real theoretical justification. 

The simplest way to arrange the condition of B = is to rely on the no-scale super- 
gravity models as in Ref. |3^]. For this, let us take the hidden sector field <j> with a Kalher 
potential 

K = -31og(</> + <f), 
and its Yukawa coupling to matter fields 

iio,i2o(0) = e~ c t Yi 26 = const. (4.6) 

This can be a consequence of an U(l) symmetry under which (f) transforms like <f> — > 
4> + ia, and then 16.16.10 and 16.16.120 are charged but 16.16.126 is not. The A terms 
associated to the Yukawa couplings are given by 

AF = -m 3/2 (0 + 0*)^Y. (4.7) 

Therefore, we obtain ^10,120 ~ m 3/2 an d B = A12Q = at the GUT scale. 

On the other hand the smallness of the B term needed for a successful soft leptogenesis 
could follow simply from a tuning among various supersymmetry breaking terms, which 
is technically natural if it is stable under sub- leading corrections. It is amusing to realize 
that the B term receives an important radiative correction due to gauge interactions of the 
right-handed (s)neutrinos. In the context of SO(W), the right-handed (s)neutrinos have a 
coupling to a heavy gauge boson X and the corresponding gaugino X which also obtains 
a supersymmetry breaking mass my 2 . Specifically, the gauge coupling N-N-X leads to 
the one-loop correction to the B term of N which is given by 

a My 
BM ps —m 1/2 M log — ^ (4.8) 

where Mx is the mass scale of the heavy gauge boson X or the B — L symmetry breaking, 
for instance. Now, putting a = 1/30, M = 10 8 GeV and M x = 10 10 GeV, we find 
B ~ W~ 2 rrii/2 which gives us the required value of B ~ 1 GeV for m 1 / 2 = 100 GeV. 
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5. Conclusions 



The motivation of our work was to understand why in many GUT models which describe 
successfully the right values of fermion masses and mixings (0, |J, |J], ||]), it is not 
possible to achieve the observed baryon asymmetry through the decay of heavy right- 
handed neutrinos. Apart from the strong hierarchy in the neutrino Yukawa couplings Y u , 
two factors are important: first linking m u and m v LR (in the simplest case they could be 
the same), which is a general 50(10) GUT relation, and then having the special feature 
of Y{[ — * 0. Starting from these conditions we reconstructed the general structure of Y v 
and the mass scales of two right-handed neutrinos which are compatible with the neutrino 
data as well as the GUT relations enforced by a certain flavor symmetry in the decoupling 
limit of the heaviest right-handed neutrinos. 

Our analysis shows that the neutrino couplings associated with two light right-handed 
neutrinos are determined by Y^- 2 ~ ~ 10~ 4 while the right-handed neutrino masses 
are of order (10 7 — 10 8 ) GeV. Conventionally, such a parameter region is far away from 
a successful leptogenesis unless a certain fine-tuning is arranged between two light right- 
handed neutrinos to resonantly enhance the resulting lepton asymmetry. However, the 
soft leptogenesis arising from the CP phases of A and B supersymmetry breaking soft 
parameters can work consistently with our picture although our parameters are in the 
strong wash-out regime of the lepton asymmetry. The basic ingredients for this to occur 
are (i) — > and (ii) a resonance condition of B ~ V. The first property protects 
the electron asymmetry which is generated by the second lightest right-handed sneutrino 
whose wash-out factor is favorably smaller. Interestingly the resonance condition requiring 
B ~ 1 GeV can be a consequence of the gauge one-loop correction involving the coupling 
of the right-handed sneutrino to the heavy GUT gaugino. For the vanishing condition of 
the tree-level B term, one may invoke no-scale supergravity, as it is difficult to achieve 
a dynamical realization of B = by arranging a specific form of the superpotential and 
Kahler potential, which requires to introduce a fine tuning of parameters in the general 
supergravity context. 
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A. Numerical form of m u 

The order of magnitudes in m u can be illustrated in the case m U3 3> m U2 > m vi , for which 
we determine numerically the approximate values of m u up to the value of t < r, Eq. ([2.4|), 
and the limit sf 3 < 0.041: 

= 0.0542 + sf 3 + 0.719t 
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y^lm = 0.102(0.600 - 0.375 s 13 ) + ^ - 0.848(0.375 + 0.600si 3 ) t 

( m ^ eT = -0.102(0.600 + 0.375 s 13 ) + + 0.848(0.375 - 0.600si 3 )t 

m„ 3 V2 

^htt = I + 0.193(0.600 - 0.376 si 3 ) 2 - (0.375 - 0.600 si 3 ) 2 t 
m V3 2 

( Tfl i 1 

v v)p!r = - - 0.193(0.600 - 0.376 si 3 )(0.600 - 0.376si 3 ) 
m VA 2 

-(0.375 + 0.600 si 3 )(0.375 - 0.600 s 13 )t 



m„ 3 2 



- - 0.193(0.600 + 0.375si 3 r + (0.375 - 0.600si 3 )^t. (A.l) 



B. Fit to a particular form of symmetric Yukawa matrices 
B.l Assumptions 



A fit of textures for up and down Yukawa matrices of the form Eq. ( p. 36 ) with Y\\ negligible 



can be made compatible with the experimental information, under the following theoretical 
assumptions: 

(a) Yukawa matrices are the only source of CP violation. 

(b) We have a supersymmetric scenario which respects at low energy the constraints of 
unitarity of the CKM matrix. 

(c) The supersymmetric corrections to the ratios ^ and ^ do not exceed the per- 
centage of error on those ratios as quoted in Table ^ 

Since we have assumed (a) and (b), we must make the fits that test the unitarity 
of the CKM matrix in the Standard Model, where all experimental information has been 



taken into account, rather than to specific experiments. We use the classical fit [39], which 
takes into account the measurements of the following flavour violating processes: 

T77T> TTTT' e fc> Am B d , Am Bs , and sin/?. (B.l) 
\Vcb\ \vts\ 

These fits are a test of the unitarity of the CKM matrix in the Standard Model. That is, 
if all the experimental inputs in Eq. ( |B.1[ ) are in agreement with the unitary of the CKM 
matrix, the statistical and systematic errors are under control on those measurements and 
there is no sensitivity to physics beyond the Standard Model in those processes then after 
the fit all these fitted quantities will agree with their input values at the 68% confidence 
level (C.L.). If some of them do not agree then there is an indication of either (i) the 
departure of the unitarity of the CKM matrix, (ii) a large correction from statistical or 
systematic errors in the experimental measurements, or (iii) a contribution from process 



beyond the standard model to the constraints Eq. ( B.l ) at the level of sensitivity at which 
the measurements and the analyses are performed. The other part of our experimental 
inputs comes from considering the following mass ratios 

m u m c m s 

, , , (B.2J 

m d m s m b 
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at the scale Mz, as well as the chiral perturbation parameter Q 

m s /m d 



Q 



a/1 - {m s /m d y 



(B.3) 



We take into account the renormalization from the low scales, at which the quark masses 
are measured or computed according to experimental data, up to Mz as 



mi{M z 



mi (mi 



for i = c, b, t ; 



Vi 



mj(Mz) 
m;(2 GeV) 



for i = u, d, s. 



At two loops in QCD they can be estimated to be r\ c = 0.56, % = 0.69, r\t = 1-06, r\ u = r] d 
= r] s = 0.65. 

The form of the Yukawa matrix Eq. ( 2.36 ) is thought as being compatible with super- 
symmetric SO(10) models for which tan/3 is large, (~ 40-50). In this case the corrections 
to quark masses are not negligible, e.g., for the b quark can be up to 20% [37|. Thus strictly 
speaking we have 



m f d 



V2M W ^ cos 0s ( 1 + e f . tan S ) , 
9 



rrij 



fu = a/2Mw/— sin/3 s (l + e fu tan/fc), 



where the parameters £f )U depend on the supersymmetric particles, such as charginos, 
neutralinos and gluinos, and yj i are eigenvalues of the Yukawa matrices. However we 
assume here that the ratios are not strongly affected by those corrections: 
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l + e c y c 



1 Vu 

Tuc Vc 



1, 
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analogously for the other mass ratios considered. 

Under the conditions of Eq. ( |2.36j ) then we expect the angles of the left diagonalization 
matrices to be given as 



'12 



s 23 



Y 22 _ Vs_ 

2^ + 1 

Vb 



m,. 



s 12 



■yd 



yd 

1 X', 



m b rsb 



N 2 ^ + 1 " 



Ud 



Vs + Vd 



111,1 



rKS 



(B.5) 



where r a {, are defined as in Eq. ( |B.4| ), and sf 2 now contains m d in the denominator, which 
is a correction to the approximate formula sf 2 = ^m d /m s . The angle is obtained, of 
course, assuming that it is small and extracted from the relation, 



s 23 



•y-d 

J 23 
•yd 
1 33 



Vs 
Vb 



\rd 

1 22 



•yd „d2 
-'33*23 



The rest of the diagonalization angles are subject to the following conditions: 
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< ^23 < S 
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If we assume that the phases of the elements Y 23 and ^23 vanish or are the same, we can 
describe then the angles defining the CKM matrix in terms of two phases, <pi, <p 2 [23]: 



si2e 



s 12 



d u i<t>i 

c 12 s 12 e 



si3e 
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(B.7) 



where then <pi and (j) 2 will be given as combination of the phases of the elements y/j, 
except for = (2,3) and for Yj\ that we are neglecting. Hence the CKM elements are 
expressed in terms of 
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ImJ = (4 3 ) 2 cf 2 sf 2 *i2 sin 0i + cf 2 sf 2 K sin((j)i 



s 12 K sin 02 



(B. 



where -ftT 



'13/ S 23- I n Eq. ( p.8[ ), the quantities that are not given by the form of the 



matrix Eq. ( |2.36| ) or by the conditions Eq. ( |B.6[) , are 

1 and <t>2 



'13' 



(B.9) 



1, we have s 23 < s i2- Also from the dependence of 



Assuming sf 3 <C s 23 and given Eq. 
Vtd /Vts ° n s i3 an d s 23 and from the fact that we are fitting to s 23 directly because we are 
effectively making V c t> = s 23 , we can fit to K and the phases 4>i and <j>2- Here we choose to 



fit also 4>i to check the level of compatibility of having <j)\ 
B.2 Method of the fit 



tt/2. 



Note that we do not know all the entries in Eq. (BJS) and we have 4 CKM parameters and 
3 mass ratios from which we can fit. Thus, instead of using a \ method, we use a Bayesian 
approach where we can obtain the combined probability distribution for K, <p 2 and <j>\, 
which is identified with the likelihood; 



C{K,(j)x,cj) 2 ) oc J Y{ f(cj\cj(K,(f>i,(/) 2 ,{xi})) X Y[ fi(xi) dxi x f (K, 0i, (f) 2 ), 



j=l,M 



i=l,jV 



(B.10) 



where f(6j\cj(K, cp 2 {xi})) is the conditional probability density function (pdf) of the 
constraints 



^3 — \Vub\) \ V c b\yus) \ Vtd\l\Vts\ and Im{J} given their dependence as functions 
ire 

{m u /m d , m c /m s , m s /m b ,Q}. 



of the texture parameters: sf 2 , s" 2 , cf 2 , sf 3 and s 23 as well as the parameters and x\ 
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We have taken the values of the constraints Cj = \V U {,\, \V c b\,V us , | V* ^ | / 1 V* s | and 
Im{J} from the most recent CKM fitter results [40] , listed in Table @. The values of the 
observables 

sin 2a, sin 2(3, sin 27, (B.H) 

that we have to compare to our fit are then the values from the same CKM fitter results. 
Here 

a = Argl-VnViJiVuV^)], (3 = Arg[-V 21 V£/(V 31 V£ 3 )] 
l = Avg[-V n V^/(V 21 Vi 3 )}, 

are the angles of the unitary triangle of the CKM matrix V . 

For our fit then the fitted values of the parameter of texture Eq. ( 2.3G| ) will be in 



agreement with the fits of the CKM fitter of the unitary triangle, assuming that the su- 
persymmetric contributions are not relevant at the sensitivity at which the parameters 
Eq. (|B.1|) are related to their SM counterpart. 



B.3 Results and comments 

In Figure [l], we show the results for the 2D probabilities of the parameter K versus (f>2 
and <f>\. Given the method used for this fit we expect at least a 95% C.L. compatibility 
with what we have fitted. If that is not the case then obviously one of our theoretical 
assumptions should be modified. However we do have a 68% C.L. compatibility of all 
of our input values with the output (fitted) information. We show in Table || the input 
values used for the fit. We have chosen to use the results of the CKM fitter collaboration 
fiofl (which provides the SM fits to the Particle Data Group). In Table ||] we show for 
comparison the different values of the unitary angles. In Table || we show the results of 
our fit. 




We can notice indeed that while the output angles and 7 are fitted in great agreement 
with the inputs, the angle a tends to be lower than the CKM fitter central value, although 
compatible at 68% C.L. This tendency could be due to the fact that the unitary triangle 
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Input Values 



Constraints 


Value ± Gaussian errors Flat errors % of error Referen. 


\v ub \ 


fo «sq+0.106\ v i n -3 
(6.b$6_ 079 ) X 1U 






[FTTH 


IT/. 1 
1 v cb\ 


( 4L61 Ia63) X 10 








IT/ 1 /It/" I 

\Vtd\/\Vts\ 


0.2003_ 0059 






u 


IT/ 1 
1 v us\ 


n nn^i r+0 00101 

u.ziz, i J.<J_n noiOO 








t r ti 

Im{ J\ 


(2.91+^f) x 10- 5 








Varied Parameters. 










m u 
m d 
m c 


0.553 ± 0.043 




7.7% 




11.3 ±2.8 




33.5% 




m s 








m 3 
m b 


0.0213 ±0.006 




28% 




Q 


22.7 ±0.8 








y d 22 


0.036 


±0.014 

















Table 2: Input values for constraints and varied parameters which are also fitted. 



Other values from the CKM fitter 



Parameter 


CKM value ±la C.L. 


±2ct C.L. 


Direct exp. value ±1<t C.L. 


±2(7 C.L. 


a 
P 
7 


22.03j£g 
59.0±|;? 


+8.0 

-17.9 

+1.69 

-1.27 

+18.0 

-7.3 


92.6.0^° 3 7 

91 OQ+1-03 
Z1 - ZO -0.99 

60.0±|| 


+27.1 

-15.7 

+2.09 

-1.96 

+62 

-39 



Table 3: Relevant information from experiments and from the CKM fitter [j40|]. The later are 
included indirectly in the fit because they are not used as constraints. 



(UT) fit itself have shown consistently during the last 5 years a difficulty in fitting the 
unitary condition itself: it = a + + 7 (which is used in such fits with respect to the direct 
measurements) as we can immediately see also in Table ^. Except for the fits of last year, 
all the angles of the UT fits have nevertheless being in agreement with that condition at 
68% C.L. But there is still the possibility that there could be a sizable beyond the SM 
contribution that could show in future analyses, therefore changing the contribution of the 
SM values to the fitted values of the unitary angles. Note also that the central value of the 
direct experimental value of a is lower than the central value of the CKM fitter and the 
errors are comparable. 

Another thing to consider, of course, is that the model based in Eqs. (2.36), ( |B.5| ) 
and ( p.7| ) would need to include corrections or modifications. Among them are the su- 
persymmetric corrections to the quark masses that should be carefully taken into account, 
or deviations from the symmetric textures. These deviations are indeed formally present 
since the symmetric structure of the mass matrices is valid just at the GUT scale and 
may get sizable modification by the RGE running to the electroweak scale. Given the 
increasing precision in the determination of the unitary triangle fits, this running should 
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Output Values 



Parameter 


Value ± errors 


K 


0.044+°"°^ 

—0.007 




1 t;qn+ 050 

— U.UoU 




-0.8777°97? 

+0.271 


viz 
viz 




\v ub \ 


(3-75^;^) x 10 3 


\ V cb\ 


(42.25I 060 ) x 10 


\Vtd\/\V ts \ 


204fi +0 ' 0096 

U.ZU4D_ 0057 


\V US \ 


997M +0-00091 

u.z^wd±_ 00081 


Im{J} 


(3.02+° -^f) x 10~ 5 


a 


(82.14i^)° 


a 


(21.83 + ?J 3 )° 


7 


«<5 = (60.3 + y } o 



Table 4: Our output values 



be taken into account. The running effects give a correction to the relation s" 2 of the form 
s i2 ~ ^/ r \/ m u/ m a where r is a parameter of 0(1) measuring the slight non-symmetry 
of the elements \Y^ 2 \ and |V 2 "|. Another possible modification is the one pointed out in 
1 24], namely allowing the contribution of YJi to become non negligible. This produces the 
same relation of s" 2 ~ l/ r \/ m u /m c with r depending on the non negligible element. 
A separate question, independent of the fit itself, is whether this fit is compatible with a 



particular realization of a horizontal symmetry, like the one proposed in [23]. 
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